Neutrino-pair emission due to electron-phonon scattering in a 

neutron star crust 



D.G. Yakovlev and A.D. Kaminker 
loffe Physical- Technical Institute 
^ ■ 194021, St.Petersburg, Russia 

o\ ■ 

, . Abstract 

Q^, Neutrino-pair bremsstrahlung radiation is considered due to electron-phonon scattering 

■ of degenerate, relativistic electrons in a lattice of spherical atomic nuclei in a neutron star 

crust. The neutrino energy generation rate is calculated taking into account exact spectrum 
of phonons, the Debye- Waller factor, and the nuclear form-factor in the density range from 
10'' g cm-3 to lO^"* g cm at arbitrary nuclear composition for body-centered-cubic and 
face-centered-cubic Coulomb crystals. The results are fitted by a unified analytic expression. 
^ I A comparison is given of the neutrino bremsstrahlung energy losses in a neutron star crust 

composed of ground state and accreted matter, in the solid and liquid phases. 



^ ■ 1 Introduction 

' The neutrino-pair bremsstrahlung radiation of electrons in the liquid and crystalline phases 

• of dense matter is one of the main cooling mechanisms of neutron star (NS) crusts. This 

radiation has been studied in a number of works (see, e.g., Itoh et al., 1989, and references 
O ' therein). Recently Haensel et al. (1966) have reconsidered the neutrino radiation produced 

-4— » ■ due to collisions of degenerate, relativistic electrons with atomic nuclei in the liquid phase, 

! e + (Z, A) — > e + {Z, A) + u + u. The authors have expressed the neutrino energy loss rate 

in a simple form through the Coulomb logarithm which depends weakly on the parameters of 
stellar matter. The Coulomb logarithm has been calculated and fitted analytically for matter 
with arbitrary nuclear composition at all densities and temperatures of practical interest. The 
\ proposed approach is invalid only in the deep layers of the NS crust (p ~ 10^^ g cm~^), where 

the atomic nuclei can form clusters and/or have nonspherical shape (Lorenz et al., 1993, Pethick 
and Ravenhall, 1995). 

In the present article, we also restrict ourselves to consideration of spherical atomic nuclei. 
We will study the neutrino-pair bremsstrahlung of electrons in the crystalline phase of dense 
matter. It is well known (e.g., Itoh et al., 1989) that this process has two channels. The 
first one is the neutrino radiation due to Bragg diffraction of electrons in a lattice ("the static 
lattice contribution" ) . The second one is the neutrino emission due to absorption or creation of 
phonons ("the phonon contribution"). It has been accepted until recently that both components 
have comparable intensities. However recently Pethick and Thorsson (1994) have shown that 
the band structure of crystals in a NS crust (the presence of an energy gap in the electron 
dispersion relation at the Brillouin zone boundaries) suppresses exponentially the static lattice 
contribution. Accordingly the phonon contribution dominates, and it will be the subject of our 
consideration. 
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The neutrino radiation due to the electron-phonon scattering has been studied earher by 
Flowers (1973) and — in the frame of the Weinberg - Salam - Glashow theory — by Itoh et al. 
(1984a). The latter authors fitted their results by complicated analytic expressions and tabulated 
the fit parameters for matter with different nuclear compositions. These results are also included 
in the new review of Itoh et al. (1996). We will show that, in some cases, the fits of Itoh et 
al. (1984a, 1996) are not too accurate. We will carry out more accurate calculations and fit the 
results by a unified formula valid for matter with arbitrary nuclear composition. In addition 
to the body-centered-cubic crystals, which have been studied traditionally in the literature, we 
will consider also the face-centered-cubic crystals. For illustration, we analyse the neutrino-pair 
bremsstrahlung in the solid and liquid phases of NS crusts composed of ground-state or accreted 
matters. 



2 Physical conditions 

Consider matter of a NS crust in the density range 10^ g cm~^3 <C p < 10^^ g cm~^. In the 
domain of 10^ g cm~^ ^ /) < /Ond — 4 — 6 x 10'''^ g cm~^, which corresponds to the outer crust, 
the matter consists of atomic nuclei immersed in a nearly ideal relativistic and degenerate gas 
of electrons. In the inner crust, where p > pND, free neutrons appear in addition to the nuclei 
and electrons. If p > lO^^g cm~^ the nuclei can be nonspherical and form clusters (Lorenz et 
al., 1993; Pethick and Ravenhall, 1995); the nuclei disappear at the bottom of the inner crust, 
at p > 1.5 X lO^'' g cm~^. 

The state of degenerate electrons in the NS crust is characterized by the Fermi momentum 
Pf or the relativistic parameter x: 

PF = ^(37^V)'/^ x = ^^m.9(^) , (1) 

rUeC KlJie J 

where fi,, is the number of baryons per electron and pi2 is density in units of 10^^ g cm~^. In 
the density range of study we have p S> 10^ g cm~^, and the electron gas is ultrarelativistic 
(x » 1). The electron degeneracy temperature is 

2 

Tf = To (Vl + a;2 - 1), Tq = « 5.930 x 10^ K, (2) 

Kb 

where fee is the Boltzmann constant. We consider the case of T ^ Tf. 

The state of nuclei (ions) is characterized by the ion-coupling parameter F = Z^e^/ {ak-^T) ~ 
0.225 X Z^/^/Tg, where Ze is the nuclear charge, a = [3/(47rni)]^/^ is the mean inter-ion distance, 
n\ is the number density of nuclei, and Tg is temperature in units of 10® K. At high enough 
temperatures, when F <C 1, the nuclei constitute the Boltzmann gas. For lower T (F > 1), the 
gas transforms into the Coulomb liquid, and the liquid solidifies into the Coulomb crystal at 
T < Tin. The melting temperature of a classical one-component plasma of ions corresponds to 
Fm = 172 (Nagara et al., 1987) and equals 



Tr 



1.32 X 10^^5/3 ( Pll) K. (3) 
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Note that the crystallization of light nuclei (H, He) is additionally suppressed by strong zero- 
point vibrations (Mochkovitch and Hansen, 1979). The most bound is the body-centered-cubic 
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(bcc) crystal (see, e.g., Brush et al., 1966). The binding of the face-centered-cubic (fee) crystal is 
weaker but very slightly. One cannot exclude that the fee crystals appear in dense stellar matter 
along with the bcc ones (e.g., Baiko and Yakovlev, 1995). We will consider the temperature 
range of crystalline crust Tu < T < Tm, where T\j ~ TpZ^^^e'^ /(hvp) is the temperature below 
which the Umklapp processes are frozen out (e.g., Raikh and Yakovlev, 1982), Tp = huj^/k^ 
7.83 X Z P12/ {A^e) K is the ion plasma temperature, ojp = \jAT^Z'^e'^n\/'mi is the ion 
plasma frequency, and m\ is the ion mass. Note that the Debye temperature of a crystal is 
Td ~ 0.45Tp (Carr, 1961). If T > Tu the electron-phonon scattering can be described using 
the approximation of almost free electrons. Then, at Z ^ 1, the main contribution into the 
scattering comes from the Umklapp processes (Yakovlev and Urpin, 1980, Raikh and Yakovlev, 
1982). 

For illustration, we will consider two models of matter in a NS crust: ground state (cold- 
catalyzed) matter and accreted matter. The effects of temperature onto the nuclear composition 
will be neglected which is justified as long as T < 5 x 10^ K (see Haensel et al., 1996). Both 
models assume that nuclei of one species are available in matter at any density (pressure). This 
leads to jumps of nuclear composition with increasing density (pressure). 

The ground-state matter in the outer NS crust will be described using the data obtained 
by Haensel and Pichon (1994) from new laboratory measurements of nuclear masses with large 
neutron excess. For the inner crust, we will use the results of Negele and Vautherin (1973) 
derived on the basis of a modified Hartree-Fock method. 

We will use the recent calculations of Haensel and Zhdunik (1990b), for describing accreted 
matter. The light elements are assumed to be burnt out in the outer layers of the NS crust while 
the nuclear composition in the deeper layers (we are interested in) is determined by a sequence 
of nuclear transmutations which occur in a sinking and compressing accretion matter and which 
depend neither on temperature nor on accretion rate. 

The comparison of the ground-state and accreted matters is given, for instance, by Haensel 
and Zhdunik (1990a, b) or Haensel et al. (1996). At p > 10^ g cm~^, the accreted matter 
is generally composed of lighter nuclei with lower Z. For instance, at p = 4 x 10^^ g cm~^ 
the accreted matter contains neutron-reach magnesium nuclei, A = 44, Z = 12 (Haensel and 
Zhdunik, 1990a, b), whereas the ground-state matter is composed of tin nuclei, A = 159, Z = 50 
(Negele and Vautherin, 1973). The above difference is important for thermodynamic and kinetic 
properties of matter. For instance, the melting temperature (^ of the accreted matter at 10^^ 
g cm~^ P ^ 10^^ g cm~'^ is about 3-10 times lower than that of the ground-state matter 
(Figure 1). 

3 General expressions 

We can use the scheme of extended Brillouin zones and the free-electron approximation for 
describing the electron states. In this case, the Fermi surface and the electron dispersion relation 
are the same as for free electrons. This approximation is violated only in the narrow regions 
of momentum space near intersections of the electron Fermi surface with the Brillouin zone 
boundaries where the dispersion relation contains the energy gap. Note that the interzone 
electron transitions which accompany the neutrino-pair emission due to the Bragg diffraction 
(Section 1) occur just in these regions. The presence of the gap suppresses these transitions 
exponentially (Pethick and Thorsson, 1994). 
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A general expression for the neutrino energy loss rate Q due to the electron-phonon scattering 
of relativistic and degenerate electrons can be obtained using the standard formalism of the 
electron-phonon interaction (Ziman, 1960). In analogy with the results of Haensel et al. (1996), 
Q can be written as 

Q = J,,9 , ^ {kBTfn,L « 3.229 x lO" pi2 Xa -r T^L ergs s^^ cm-^, (4) 

where Gp = 1.436 x 10~^^ ergs cm'^ is the Fermi weak-coupling constant, Xa is the mass 
fraction contained in nuclei, L is a dimensionless function of the parameters of stellar matter 
which has meaning of the Coulomb logarithm in the liquid phase. Furthermore, C\ = Cy + C\ + 
2{C'y + C'a), where Cy and Ca are the vector and axial-vector constants of weak interaction, 
respectively. We have Cy = 2sin^0^y + 0.5 and Ca = 0.5, for the emission of the electron 
neutrinos (charged + neutral currents); and C'y = 2sin^0vK — 0.5 and C'a = —0.5, for the 
emission of the muonic and tauonic neutrinos (neutral current only). In this case, 9^ is the 
Weinberg angle, sin^ 6w — 0.23. Numerical estimate of Q is obtained taking into account the 
contribution of the electron, muon and tau neutrinos which yields « 1.675. 
The general expression for L in Equation (^) can be presented in the form 




y2 



(5) 



where 



— = 0.0128Tq — , z = — 



w + z w 
+ 



(6) 



.(1 - e-^)(e"'+^ - 1) (e^ - l)(e^-^ - 1 

In this case, s= 1, 2, 3 enumerates three phonon modes in a Coulomb crystal. The integration 
in Equation @ is carried out over possible momenta of electrons p and p' on the Fermi surface 
before and after scattering; S = 47rpp is the Fermi surface area. The quantity hq = p — p' 
is an electron momentum transfer in a scattering event and y = hq/(2pp). Furthermore, k, 
uJs{k) and es{k) are, respectively, the wave vector, frequency and polarization unit vector of a 
phonon excited or absorbed by an electron: itfc = q — K, where K is an inverse lattice vector 
chosen in such a way that k belongs to the first Brillouin zone; the neutrino-pair momentum 
can be neglected here. The case of k = q {K=0) corresponds to a so called normal scattering 
process, while the case of / q {K ^ 0) corresponds to an Umklapp process. The quantity 
e(g) in Equation (^) is the longitudinal static dielectric function of the electron gas: it takes 
into account the electron screening of the Coulomb potential of a nucleus; /(g) is the nuclear 
formfactor to allow for a finite nuclear size. Two terms under the integral in ii{z^ describe 
the neutrino-pair radiation due to emission and absorption of phonons. The dimensionless 
integration variable is w = £y/{k^T), where Ei, is the total energy of a neutrino-pair. Finally, 
W{q) is the Debye-Waller factor. Baiko and Yakovlev (1995) fitted this factor for a Coulomb 
crystal (with a mean error of about 1%) by the expression 2W{q) = ay"^, where 



o = .0 (i«-.e-"»«' + <„_.) , = » 1.683;^. (7) 
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Here, u„ = ((cj/wp)") is a frequency moment of the phonon spectrum. Brackets mean averaging 
over phonon frequencies and polarizations: 

ifsik)) = f dfe/,(fc), (8) 

Vb is the Brihouin zone volume. One has u^^l'^^ = 2.800, 7x^2^^^ = 12.998 for the bcc lattice 

(e.g., Mochkovitch and Hansen, 1979), and u^Ii'^^ = 4.03, ^^2^^ = 28.8 for the fee one (e.g., 
Baiko and Yakovlev, 1995). According to Equation (|5[), the Debye- Waller factor suppresses the 
electron-phonon interaction at large momentum transfers q, and weakens the neutrino emission. 
The factor is important at large ion vibrations in the lattice — near the melting point (where 
the thermal vibrations are especially strong) and at high densities (where zero-point vibrations 
are huge). Under the conditions of study (Section 2) this factor is most significant. 

In calculations based on Equation (||) we will use the static dielectric function of the degener- 
ate relativistic electron gas e{q) derived by Jancovici (1962). In our case, the electron screening 
is weak and has practically no effect on the results. Nevertheless we will take it into account. 

We will use the nuclear formfactor appropriate for spherical nuclei with a uniform proton 
core: 

3 

fig) = [sin(grc) - qrc cos(grc)] , (9) 

where, Tc is the core radius (which can be noticeably smaller than the total nucleus radius if 
the nucleus is neutron-rich). The formfactor suppresses the electron-phonon scattering at large 
momentum transfers q. One has f{q) = 1, for point-like nuclei. 

It is convenient to fit H(z) by an analytic function. One can easily see that H{0) = 1, and 
H{z) = 21 e~^/(1607r^) at 2: ^ 00. We have calculated H{z) numerically for intermediate 
values of z. The calculated values and the asymptotes are reproduced (with the mean error of 
1.7 %) by the fit 

H{z) = (0.2968 + V0.01633z2 + 0.70322)^ exp (5.5032 - Vz"^ + 5.50322) _ (^q) 



4 Calculations and fitting of the results 

Evaluation of L reduces to a 4-fold integration in Equation (|5|) over positions of electron momenta 
on the Fermi surface before and after scattering. For some parameters of stellar matter, the 
integration has been carried out numerically by the Monte Carlo method. We have used the 
code developed by Baiko and Yakovlev (1995) for calculating analogous integrals which determine 
the thermal and electric conductivities of electrons due to electron-phonon scattering. We have 
modified the code for calculating L. Typical calculation error for 10^ Monte Carlo configurations 
(choices of positions of electron momenta on the Fermi surface before and after scattering) has 
been 6 % (at the 3a level). 

In order to simplify calculations and analyse the results, we have used also the simplified 
method of integration of the equations similar to (^), which was used earlier by Yakovlev and 
Urpin (1980), Raikh and Yakovlev (1982), Itoh et al. (1984b, 1993), and Baiko and Yakovlev 
(1995) in calculations of the thermal and electric conductivities. The method is based on the 
fact that the main contribution into the electron-phonon scattering comes from the Umklapp 
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processes. The Fermi surface is intersected by many boundaries of the Brillouin zones. Then, 
for the Umklapp processes, in Equation (|5|) one can approximately set 

Y^[qes{k)]^H{z)z'^ ^ qH^G{t), (11) 



G{t)=t-^{H{z)z-^), (12) 
where brackets denote the averaging (^). As a result, from Equation (|5|) we obtain 

L= aotG{t)P. (13) 

Here, 



,,2 



yo 



y = hq/{2pp), and the lower integration limit yo = ^Qo/C^Pf) = (4Z)^^/^ is chosen in such a 
way to include the Umklapp processes but exclude the normal ones. The minimum momentum 
qo = (Gvr^ni)^/^ is set equal to the radius of the sphere whose volume equals the Brillouin zone 
volume. 



Calculation of L from the approximate expression (13) reduces to determination of the two 
functions, G and P, that is much simpler than the exact Monte Carlo integration. Thus the 
calculations have been done mainly by the simplified method whose accuracy has been controlled 
by the Monte Carlo runs at some parameter values. In all the cases, the results have coincided, 
within the Monte Carlo error bars. 

The function G depends on the dimensionless temperature t, given by Equation (^), and on 
the lattice type. If t ^ 1, from Equations (^) and (|l2|) we obtain the asymptote G{t) = U-2, 
where U-2 is a frequency moment of the phonon spectrum (Section 3). When t ^ we have 
G{t) ~ t. We have calculated G{t) from Equation (|l^ for the bcc and fee lattices using 
the approximate but highly accurate method proposed by Mochkovitch and Hansen (1979) for 
averaging the expressions similar to (^) over the Brillouin zone. The results can be fitted (with 
about 1.5 % mean error) by the function 



where af''''^ = 0.06423, 6^''^ = 0.1151, 6^^''^ = 5.92 x 10"^ af,^"^^ = 0.001479, fej^'^^ = 0.2381, 
^2^'^'^^ = 5.161 X 10~^. This expression reproduces also the above asymptotes. 

The function P is given by a simple one-dimensional integral, which is insensitive to the 
crystal type but depends on three parameters: atomic number Z (through the lower integration 
limit, yo)) the Debye-Waller parameter a = a{t) (see Equation (|7|)), and the parameter r] = Vc/a 
which accounts for the finite size of atomic nuclei in the nuclear formfactor (P) (a is given by 



Equation (|)). We have evaluated P from Equation (yD at Z = 20,40,60, a = 0.04, 0.12, 0.4, 
1.2, 4, 12 and rj = 0, 0.10, 0.2, 0.3, 0.4. Analysing the parameters of matter in a NS crust at 
10^ g cm~^ ^ P ^ 10^^ g cm~^ (Section 2), one can show that the selected values cover all the 
domain of possible parameter values. The results are fitted by the expression 

-N 

(16) 



P(Z,a,r?) = ^ 



1 I {l8nZr/' Fi2 
5N Fn ' 
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where 2pYrc/fi 



(187rZ)i/377, N = 3.781 - 0.02322 Z, and the functions Fq and Fi are given by 



a 

c1 



-0.3462Q 



0.6449 



l + a 

:-[i„W-o.51 e-o + i^ 

2 (1 + aj'^ 

1 



In 



Fi = ^ [(l + n)e-" -(l + a) e""] - 0.3949 



1 + n 
l + a 
e" 



- 0.4230 

0.4726O 



a 



(1.5 + a)3/4 



[hi(s) - 5.213] e" 



a (2 + 4n + 
(1.5 + a)4 



In 



1 +u 
1.5 + a 



0.1 



(17) 



(18) 



Here, s = (4Z)~^/^ + 0.00563, u = sa. The mean error of the fit (16) is 1.3 %, and the maximum 
error of 3.9 % takes place at t=0.12, Z=60, r/=0.3. 

Note a simple asymptote L ~ 0.5 which is valid in the limiting case when t ^ 1, a ^ 1 and 
a^Q <C 1. However this asymptote cannot be actually realized in NS matter for the conditions 
of study (Section 2). In the opposite case, when t <C 1, we obtain L = aot^n_2-P/ao, where P 
should be calculated at a = 0.5aoti-i (see Equation (|^). 



5 Discussion 



Equations (§) , (|^) , (ITsI) and (|T6|) allow us to evaluate easily the neutrino energy loss rate Q due 
to electron-phonon scattering for any model of stellar matter in the parameter range of study 
(Section 2). 

Figures 2-4 display L as a function of the dimensionless temperature t. The maximum 
values of t correspond to the melting temperature (^). Following Itoh et al. (1984a), the proton 
core radius rc in the nuclear formfactor (|^ has been taken as: Tc = 1.15A^/^ fm at p < 4 x 10^^ 
g cm~'^; and rc = 1.83Z^/^ fm at higher p. If p < 10^^ g cm~'^ the nuclear size effect is 
unimportant and one can set r/ = 0. 

Figure 2 corresponds to the ground-state matter (Section 2) at p = 4x 10^^ g cm~'^ (A = 159, 
Z = 50, Negele and Vautherin, 1973). The nuclear size effect is seen to decrease L and Q by (20- 
30) %. According to Equations (^) and (0), the nuclear formfactor decreases the integrand at 
large ?/ ~ 1 (i.e., at large-angle scattering), and the effect is more pronounced with increasing rj. 
However the integrand is already small at these y due to the Debye- Waller factor (which reduces 
the scattering with large momentum transfers) and due to the factor [l+(2y^/(l— y^)) ln(y)]. The 
latter factor describes suppression of backward scattering of relativistic electrons (see Haensel 
et al., 1996) and vanishes at y = 1. This explains weak dependence of L on rj. 

The temperature dependence of L in the bcc and fee crystals is similar. The L{t) curve for 
the bcc crystal has formally the maximum, just as the bcc curve, but this happens at higher t, 
above the melting temperature. The maximum of L(t) for the fee crystal is shifted to lower t 
(and seen in Figure 2) owing to a softer phonon spectrum (see Baiko and Yakovlev, 1995). 

Let us emphasize that the neutrino emission due to the electron-phonon scattering in the 
bcc crystal was calculated earlier by Itoh et al. (1984a). These results are also included in the 
new review article of Itoh et al. (1996). The authors calculated the function -^phonon 

related to 

our function L as -Fphonon = 2L/3. The calculations of Fphonon were carried out for the nuclei of 
^He, 12c, 1^0, 20Ne, ^^Mg, 2831, 323^ 40ca, 56pg 10^ g cm-^ < p < 10^2 g cm-^ and also for 
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the ground-state matter at 10^ g cm~'^ ^ P ^ 10^^ g cm~^. The results were fitted (Itoh et al. 
1984a, 1996) by very comphcated expressions which require large tables of fit parameters. Our 
fits (Section 4) are much simpler and valid for any nuclear composition. For comparison, Figure 
2 presents the temperature dependence of L{t) (long dashes) for the bcc crystal derived using 
the fits of Itoh et al. (1984a, 1996). The results are somewhat different from those obtained in 
the present article. The difference (which is sometimes rather large) takes place also for other 
values of /9, T, and for other elements (Figures 3 and 4). 

This difference is likely to come from not too accurate fits of Itoh et al. (1984a, 1996). In order 
to confirm this statement let us compare the equations of Itoh et al. (1984a) and the equations 
of the present article. Itoh et al. (1984a) used another method of calculation taken from the 
earlier article of Flowers (1973). For comparison, we express the integral ( |l4|) in Equation ( |l3|) 
as a three-dimensional integral over electron momentum transfers dq. Flowers (1973) and Itoh 
et al. (1984a) replaced the latter integration by summation over possible inverse lattice cells 
where q could appear. The integration over any cell of the sum was done by replacing q = K. 
In addition, just as in the present article, the normal scattering processes were neglected, i.e., 
the contribution from the central cell was excluded. For comparison, we have also calculated L 
by the same method of summation over inverse lattice cells. The results (short-dashed curves in 
Figures 3 and 4) are in good agreement with our original results (for chosen nuclei and densities) 
but differ from those given by the fits of Itoh et al. (1984a, 1986) by a factor of several. 

Figure 5 demonstrates the neutrino energy loss rate E = Q/p per unit mass (ergs g~^) as 
a function of density at several values of temperature for the ground-state and accretion matters 
(Section 2). Note that the selfconsistent models of accreted matter correspond to temperatures 
T ~ 10^ K (Miralda-Escude et al., 1990). We have extended the accretion model to higher T 
for illustrative purpose, to demonstrate the effect of variation of nuclear composition onto the 
neutrino emission rate. For the displayed parameters, the matter can be either in the liquid or 
in the solid phase. The emissivity E in the solid phase has been calculated using the expressions 
of Section 4 (for bcc crystals). In the liquid phase, we have used the formulae of Haensel et 
al. (1996). Nuclear composition varies at certain values of density as described in Section 2 
(Figure 1). In addition, there occur liquid-solid phase transitions at certain p which are easily 
deduced from Figure 1. For all T displayed in Figure 5, the matter of sufficiently high density 
becomes solid. With decreasing T, the crystallization boundary shifts towards lower densities. 
The accreted matter solidifies at somewhat higher p since it contains nuclei with lower Z. In 
addition, if T = 8 x 10^ K and p increases in the accreted matter and becomes higher than 
p > 10^^ g cm^^ there occurs a series of transitions from liquid into solid and back. This effect 
is associated with strongly nonmonotonic melting temperature due to rapid changes of Z (Figure 

!)• 

According to Equation (Q) , the dependence of the neutrino energy loss rate E on p and T is 
determined by the factor (XaZ'^/A) r| L. If XaZ'^/A and L were fixed, the emissivity E would 
be displayed by horizontal lines in Figure 5. Actually, however, E varies with p. The curves 
are seen to suffer jumps associated with changes of the nuclear composition and crystallization. 
The jumps due to variations of the nuclear composition are much weaker. In a crystallization 
point, the function L and the emissivity E decrease steeply by a factor of 1.5-3 but the general 
character of the density dependence of E in the liquid phase and just after the crystallization is 
the same. This is because the neutrino-pair bremsstrahlung by electrons due to the Coulomb 
scattering in a liquid is similar to that due to the high-temperature electron-phonon scattering 
(T > Tp). Indeed, the function L is rather smooth either in the liquid (Haensel et al., 1996), or 
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in the high-temperature crystal, and it behaves as cx {XaZ^/A)T^, in the first approximation. 
The crystalhzation jump of E can be somewhat smoothed if we additionally take into account 
the exponentially suppressed static lattice contribution (Pethick and Thorsson, 1994). However, 
our preliminary estimates based on the approximate expressions given by Pethick and Thorsson 
(1994) and Itoh et al. (1996) show that the smoothing is not significant. 

With increasing p inside the crystalline phase (at a fixed temperature T), the ratio t = T/Tp 
generally decreases, and the low— temperature phonon scattering becomes important. In this 
regime, according to the results of Section 4, the function L = aQf?u^2P/O'0, and P depends 
weakly on the parameters of stellar matter. Then one has a strong (~ t'^) suppression of 
the neutrino generation. The generation rate becomes a sharper function of T and depends 
explicitly on p. In the first approximation, E oc T^(yl/Z)^/^ -^a^ p~^l^ . The suppression leads 
to a noticeable decrease of E with p (Figure 5). The dependence of E on nuclear composition 
becomes much weaker, than in the liquid or high-temperature crystal. This explains increasing 
similarity of the emissivities -E'(p) for the ground-state and accreted matters in Figure 5 with 
decreasing T and/or increasing p. At T < 5 x 10^ K the difference of the emissivities becomes 
negligibly small. 

As follows from Figure 5, the neutrino emissivity E is much more sensitive to the state of 
matter (liquid or crystal) than to nuclear composition (accreted or ground-state matter). For 
instance, it is easy to calculate the emissivity E for matter composed of iron nuclei ^^Fe and 
show that it differs slightly from the curves presented in Figure 5. 

6 Conclusions 

We have calculated (Section 3) the neutrino energy loss rate due to electron-phonon scattering in 
the Coulomb bcc and fee crystals of spherical atomic nuclei at densities 10^ g cm~^ ^ P ^ 10^^ 
g cm""^ and temperatures Tu < T < Tm in matter of arbitrary nuclear composition (Section 2). 
The results are fitted by simple expressions (Section 4). Combined with the results of Haensel et 
al. (1996), which describe the neutrino emission due to the bremsstrahlung in the liquid phase 
(Tin < T < Tp ) , the equations obtained allow one to evaluate easily the neutrino energy loss rate 
due to the bremsstrahlung process in the wide temperature range Tu ^ T < Tp, most important 
for applications. The main properties of the neutrino energy losses are analysed in Section 5 for 
two models of matter in a NS core (described in Section 2): ground-state and accreted matter. 
Note that at very low temperature (T <C Tm) the neutrino emission due to the electron-phonon 
scattering may appear so weak that another neutrino generation mechanism becomes dominant 
— the Coulomb scattering of electrons by charged impurities. The corresponding neutrino 
energy loss rate can be easily calculated, for instance, from the results of Haensel et al. (1996). 
However the neutrino emission at these low temperatures is expected to be too weak to be of 
practical interest. 

The results of the present article can be useful for numerical modelling of various processes 
associated with thermal evolution of NSs. First of all, we mean cooling of young NSs (of age 
< (1 - 10^) yrs), where thermal relaxation is not yet achieved (Lattimer et al., 1994). The 
relaxation is accompanied by a cooling wave which goes from the stellar interior to the surface. 
Corresponding variations of the surface temperature are, in principle, observable. The dynamics 
of thermal relaxation is quite sensitive to the properties of matter in the NS crust, especially to 
nuclear composition and neutrino energy losses. 
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Figure 1: Melting temperatures of the ground-state (solid line) and accreted (dash line) matter. 
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Figure 2: Quantity L versus t = T/Tp for the bcc (lower curves) and fee (upper curves) crystals 
composed of ^^^Sn nuclei at density p = 4 x 10^^ g cm~^ with account for the finite nuclear size 
(77 = 0.15, solid lines) and for point-like nuclei (ry = 0, short dashes). Long dashes — fits of Itoh 
et al. (1984a). 
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Figure 3: Quantity L versus t for the bcc crystal composed of °"Fe nuclei at p = 10 g cm~ . 
Solid line — present calculation, long dashes — fit of Itoh et al. (1984a, 1996), short dashes — 
present calculation using the method of Itoh et al. (1984a). 
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Figure 5: Neutrino energy loss rate due to the bremsstrahlung produced by electron-nucleus 
(r > Tin) or electron-phonon scattering (T < Tm) versus density for the ground-state (solid 
lines) and accreted (dash lines) matters at several values of T (values of Tg are shown near the 
curves) . 
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